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Abstract 



, We investigate the exact results of the Navier-Stokes equations using the methods 



developed by Polyakov. It is shown that when the velocity field and the density are 
Oh! 

(— I \ not independent, the Burgers equation is obtained leading to exact N-point generating 

•rH . functions of velocity field. Our results show that, the operator product expansion has 

X 

' to be generelized both in the absence and the presence of pressure. We find a method to 

determine the extra terms in the operator product expansion and derive its coefficients 
and find the first correction to probablity distribuation function. In the general case 
and for small pressure, we solve the problem perturbatively and find the probablity 
distribuation function for the Navier-Stokes equation in the mean field approximation. 
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1- Introduction 

A theoretical understanding of turbulence has eluded physicists for a long time. A sta- 
tistical theory of turbulence has been put forward by Kolmogorov [1], and further developed 
by others [2-4] . This is to model turbulence using stochastic partial differential equations. In 
this direction, Polyakov [5] has recently offered a field theoretic method for deriving the prob- 
ability distribution or density of states in (l-l-l)-dimensional in the problem of randomaly 
driven Burgers equation. The importance of the Burgers equation is that, it is the simplest 
equation that resembles the analytic structure of the N-S equation, at least formally, within 
the scope of applicability of the Kolmogorov's arguments. 

Polyakov formulates a new method for analyzing the inertial range correlation functions 
based on the two important ingredients in field theory and statistical physics namely the op- 
erator product expansion (OPE) and anomalies. He argues that in the limit of high Reynold's 
number because of existence of singularities at coinciding point, dissipation remains finite 
and all subheading terms give vanishing contributions in the inertial range. By using the 
OPE one finds the leading singularities and can show that this approach is self-consistent. 

Here we generlize Polyakov's approach to the Burgers turbulence [5] in the presence of 
the pressure (i.e. the Navier-Stokes equation). We consider two distinct situations, firstly 
when the velocity field and the density are not independent, and secondly when the gradient 
of the pressure is small. In the first case we show that by a change of variables one can 
transform this problem to the Burgers equation and fllowing [6] we find the exact N-point 
generating function of the velocities. For the second case we calculate the effect of pressure 
perturbatively and find the assymptotic behavior of probablity disturbuation function in the 
first approximation. 
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On the other hand, using the results of [6] , we calculate the OPE coefficiente proposed by 
Polyakov, and show that we have to correct the OPE both in the presence and the absence 
of the pressure. We write the generlized OPE and find the first correction to the PDF which 
was initially found by Polyakov. 

The paper is organised as follows. In section two we show that when the velocity and 
density are not independent the governing equation reduces to the Burgers equation. In 
section three we find the N-point generating functions for the velocity field and show that 
we have to generhze the OPE proposed by Polyakov. We find the OPE coefficients and 
derive the first correction to probablity disturbuation function. In section four we consider 
the effect of pressur perturbatively and find the probablity disturbuation function in the 
mean field approximation. 

2- The Compressible fluid and the Burgers equation 

We begin our discussion with the system describing compressible flows in an ideal (poly- 
tropic) gas ignoring dissipative effect [7], viz. 

Pt + vp^ + fW^^Q (1) 

p{vt + vv^) + ^ Q (2) 
p — kp"^, s — const. (3) 

Here p, v and p are the density, velocity fleld and pressure, respectively. The s is entropy, 

c 

assumed to be constant, 7 = the ratio of specific heats and /c is a constant. 

Also in inertial range, we ignore the driving force, which pumps energy in the large scale 
of system. The system of eqs.(l-3) has essentialy two variables (p and v), it is nonlinear and 
difficult to handel in complete generality. Nevertheless there has been considerable analytic 
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interest in this system. In particular, some progress can be made by seeking simple wave 
solutions where one of the dependent variables is a functions of the others. Rewriting eq.(l- 
3) in terms of v and p only by introducing the square of the speed of sound, — (l^)^^^^ = 
kjp'^'^, we have: 



Pt + VP:c + pVa; = (4) 
0?' { ) 

Vt + vv^ H p^ = (5) 

P 

Now, we assume that v — v{p),so that eqs.(4) and (5) changes to the following equations: 

pt + {v + pvf)px = (6) 

Pt + {v + — )p, = (7) 

Here the prime denotes differentiation with respect to p. This system of linear algebraic 
equations in pt and px has a non-trivial solution provided the determinant of the coefficient 
matrix vanishes so that: 

vf^±-^±ao{-^)'^- (8) 
P Po P 

The system (6-7), then reduces to one of the equations: 

Pt + {v± a)px = (9) 

where 

v{p) = r ^dp = -^{a{p) - ao} (10) 

We restrict our attention to waves moving to the right thus choosing the plus sign in eqs. (8) 
and (9). The corresponding equation for v follows easily from multiplication of eq.(6) by 



vf{p) and writing the result in terms of v via eq.(lO). Therefore we obtain: 

7+1 

vt + (oo H ^v)v^ = (11) 

Cole generlizes the eq.(ll) in the presence of small kinematical viscosity and have found the 
following equation for v [8]: 

7 + 1 

Vt + (oo H ^v)v^ = uv^^ (12) 

This equation can be transformed into the Burgers equation by simple change of variables. 
We define u — qq + and in the presence of some driving force in large scale eq.(12) 
transfomed to: 

Ut + uu^ = vu^x + /(x, t) (13) 



where f{x,t) is a Gaussian random force with the following correlation: 

< fix, t)f{x,t) >^ {^^fk{x - x)S{t - t) (14) 



The transformation, u{x,t) — —Xdxh{x,t) maps eq.(13) to the well knwon Kardar-Parisi- 
Zhang equation [9] , 

dth = ud,,h + ^[d^h]^ + Mx, t) (15) 

which gives the general features of many complex system [9,10]. It is noted that the parame- 
ter A that appears in the above transformation is not renormalized under any renormalization 
procedure [11]. 
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3- The Exact N-Point Generating Functions and the Generelized OPE 

To statistical description of eq.(13) following Polyakov [5] consider the following gener- 
ating functional 

N 

-^iv(Ai, A2, . . . Xn, xi,... xn) =< exp(^ Xju{xj, t)) > (16) 

Noting that the random force f{x, t) has a Gaussian distribution, Zn satisfies a closed 
differential equation provided that the viscosity u tends to zero: 

Zn + Y.^j-§^.(Y.-^) ^^Hxi - Xj)XiXjZN + Dn (17) 

where is : 

Dn ^ i^'^Xj < u"{xj, t) exp ^ Xku{xk, t) > (18) 

and 

k{xi - Xj) = {^^-^fk{xi - Xj) (19) 

To remain in the inertial range we must, however, keep p infinitesimal but non-zero. 
Polyakov argues that the anomaly mechanism implies that infinitesimal viscosity produces a 
finite effect. To compute this effect Polyakov makes the F-conjecture, which is the existence 
of an operator product expansion or the fusion rules. The fusion rule is the statement 
concerning the behaviour of correlation functions, when some subset of points are put close 
together. 

Let us use the following notation; 

Z{Xi, X2, . . . ,xi, . . . Xn) =< exi (xi) . . . cajv {xn) > (20) 

then Polyakov's F-conjecture is that in this case the OPE has the following form, 

d 

eAi(a; + y/2)eA2(a;-y/2) = A{Xi, X2,y)ex,+X:,{x) + B{Xi, X2,y)-^ex,+\:, + o{y^) (21) 

6 



This implies that Z^^ fuses into functions ^at-i as we fuse a couple of points together. 
The F-conjecture allows us to evaluate the following anomaly operator (i.e. the D^r-term in 
eq.(5)), 

a\{x^ — \rc&v{\u (x) exp(AM(a;)) (22) 



which can be written as: 
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As discussed in [3] the possible Galilean invarant expression is: 

axix) = a{X)exix) + p{X)^ex{x) (24) 

Therefore in steady state the master equation takes the following form, 

d d - 

^(— -/3(Aj))— ZAr-^A;(xi-Xj)AiAj-Zjv = 5^a(Aj)Zjv 



' dXj dxj 



/3(A) = m + j (25) 



Let us consider following correlation for / in /c-space as follows [12]: 

< f(k, t)f(k',t) >= ^ko{^^rS{k' - ^)5{k + l^)5{t - t) (26) 
therefore we obtain: 

</>,t)/V,^') >^hC^?cos{'^)5{t-t) (27) 
In the inertial range where xi — Xj « L, we can expand the r.h.s. of eq.(27) and find 

K{x, - = ^0(^)^(1 - ^^^^) (28) 
Therefore we find the following explicit form of Z2 for k{xi — xj) given by eq.(28): 

Z^iiiy) = e'^^^y^'" (29) 



and the following experssion for density of states as the Laplace transform of ; 

/•c+ic» All 

W{u,y)^ ^.e-^-ZM (30) 

where 

II ^ 2{Xi - X2), y^xi-X2. 

It can be easly shown that with the following definition of variables Polyakov's master 
equation (i.e. eq.(25) with the scaling conjecture[5] is: 

^^-^ + - + ^3/^3)^/3 = (31) 

where 

/a = (A1A2A3) ^'^Z^ 
xi + X2 + X3 X2 + 

Vl = ^ ^2=2^1 ^ yz^X2- X2, 

A1 + A2 + A3 3 A2 + A3 
/^i = ^ /^2 = -(Ai — ) /X3 = 2(A2-A3) (32) 

Now we set as; 

f-i = ^t' 9-i{^2y2, /^sl/s) (33) 

inserting in eq.(15) results in 

^3(/^2Z/2, /^3, 1/3) = e'^^^-^^+^^y^^"" (34) 

and 

^2 = ^3 = -5/4 
Now if we use the following transformation; 

Xi^- X2+ X^^- . . .Xn 



yi 



N 



= — — 

+ a;4 + . . . Xiv 

ya - -2 

and un — xn-i — xn 



and 



and 



Hi = 






N 


1^2 = 


N-1 




N -1 


1^3 = 


N-2 


I^N = 


2{Xn- 



Ai + As + . . . + A 



N 



N 

A2 + A3 + . . . + Aat, 

N^l ] 

r . A3 + A4 + . . . + AaT , 

N^2 ^ 

I — Ajv) 



we get the following partial differential equation for /jv; 

52 Q2 7 + 1 

^^-^ + ■■■+ a a }/iV - (^l^)'(y2/^2 + . . . + VNPiNffN = 

which is solved by: 

The parameter for the iV-point generating functions can be determined by the 
quirment that Zjv must be finite in the limit Aj — > 0. Therefore we find: 

2N-1 



2N 

It follows that the behavior of N-point correlation functions of v is: 

N 



k=0 



where a^.^^ are some constants. 
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Now we try to find the OPE coefficients in eq.(21). 

Tliis can be done by using tlie Z3 and put X2 and 2:3 close togetlier. Here we calculate 
the OPE coefficiens for the case 7 = 1 i.e. the ordinary Burgers turbulence and for arbitrary 
7 we need only small modifiation. and Z2 are given by eqs.(32) and (29). In eq.(32) we 
take = X2 — 2e and it is easy to show that in the limit e — > we find: 

Z, ^ iX,X2Xsf%i.2l^sr/'e'^'''^^^^^^^ (41) 

where A3 = — (Ai + A2). Now we can expand the exponent of eq.(41) in termes of e, therefore 
we find: 

Zs = (AiA2A3)^/^(//2//3)-'/V/3('^^^^+'^32/3)^''(l-^(2Ai-5A2-5A3) 
( 2(Ai - A2 - A3))'/'(a;i - X2f^^{6Xi + I3A2 - I9A3) 
( 2(Ai-A2-A3))^/'(x-i-a;2)^/'e + 0(e)^) (42) 

Comparing with eq.(21) we find that: 

A{Xi, X2, e) = (-AiA2(Ai + A2))'/'(^Ai(Ai + 2A2))-t (43) 

^(A.A2,e) = -^i^(^^^ii^ (44) 

C(Ai,A2,e) = ^A(Ai,A2,e) (45) 

Where C(Ai, A2, e) can treated as the third terms in the OPE of eq.(21) and are given in the 
following modified OPE: 

d 

exi{x + e/2)ex^{x- e/2) = A{Xi,X2,e)ex^+\^{x) + B{Xi,X2,e)—ex^+\^ 

+ C(Ai,A2,e)|^e,,+A,+0(e2) (46) 
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Now one can show that the modified OPE leads to following equations for ax{x) and Z2: 

d d 

ax{x) = a{X)ex{x) + P{X)—ex{x) + 7(A) ^e^lx) + • • • (47) 

(d^ - -)dyZ, + cfid^Z^ - fJ,VZ2 = (48) 
A* 

where we have used /9(A) = ^ , a — according [4], and scaling arguments show that 7(A) 
has the following form: 

7(A) = cA (49) 
The eq.(49) have the following assymptotic behavior: 

Z^d^y) = e'Ti('^2')'^'-§('^f) (50) 

It appers that in the limit c ^ 0, we find Polyakov's result for Z2. Indeed our results shows 
that we can find the exact Z2 for the Burgers equation by calculating the exact form of OPE, 
proposed by Polyakov by methods disscussed above. 

4- Perturbative Calculations of the Pressure 

To consider the general situation we have the following set of equations describing the 
compressible fiuid: 

Pt + vdxP + pdxV = (51) 

vt + vd^v + ^ = udlv + f(x, t) (52) 
P 

p = kp'^ (53) 

By using the eq.(53) we find: 

Pt + vd^p + pd^v = (54) 
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vt + vd^v + k-tp'-''d,p = vdlv + /(x, (55) 

where the v is the viscosity. We rewite the eqs.(54) and (55) for different points Xi and 
X2 and multiplying the resulting equations by p{x2)i\ip{x2)i p{xi) and A2p(a;i) respectively, 
and all of them to 

and averaging with respect to the external force we find: 

E(^-/9(A,-))^^2-EM^i-^.-)A.A,^2 + F2 = ^a(A,-)^2 + A2^2 

j=l ^^j '^^3 j=l j=l 

(3{X) = m + \ (57) 

where 

Z(Ai, A2, xr, X2) =< p(xi)p(x2)e^^^(^i)+^^^(^^) > (58) 
F, =< A;7p(xi)p(:r2)e^^-("^)+^^"("^)[Aip^-2(xi)9,,p(xi) + X2P''-\x2)d,,p{x2)] > (59) 

and 

2 ~ 1 

A2 = meaniY, P{^i)^—:dxiP{^i)) (60) 
For small pressure we can approximate the eq. (59) in the mean field approximation as: 

F ^ -TjZ, T] ^mean{-k-f[Xip^~'^{xi)d^^p{xi) + X2p''~'^{x2)d^2pix2)]) (61) 

By using /3(A) = j and a (A) = we find the following asmyptotic expression for Z2: 

Z2 = ei('^f)'^'-(^^+^)('^J')"'^' py^oo (62) 
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